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We review recent developments in real-time simulations of SU(2)-Higgs theory near the electroweak transition 
and related topics. 



1. Nonequilibrium phenomena 

Nonequilibrium quantum field theory is impor- 
tant in various physical situations, such as do- 
main formation during cosmological phase tran- 
sitions and the properties of the electroweak tran- 
sition in the early universe, or the QCD transi- 
tion in heavy ion collisions. We concentrate in 
this talk on the SU(2)-Higgs model and the elec- 
troweak transition, which is relevant to theories 
of baryogenesis JlJ . 

The description of nonequilibrium phenomena 
involves real time, as opposed to the imaginary 
time which is so useful for Monte Carlo computa- 
tions in quantum field theory. Consider the time 
dependence of an observable O, 

(0(t)) = TT P e lHt Oe- lHt . 

We know how to turn this operator expression 
into path integral form. It leads to phase factors 
exp(±ii9), which are hopeless for Monte Carlo. 
Very hard is also the anaytical continuation of an 
imaginary time expression back to real time when 
the data for (0(—ir)) have errors, especially at 
'large' times. A possible way out of these difficul- 
ties is the classical approximation. 

2. Classical SU(2)-Higgs model on a spa- 
tial lattice 

The action S — J dt (K—W), is given schemat- 
ically by the kinetic energy 
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Here time and Aq are in physical units, but ev- 
erything else is in lattice units a = 1, with the 
exception of the explicitly indicated lattice spac- 
ing a. 

Starting from an arbitrary initial field config- 
uration the system will evolve in time over a 
region in phase space compatible with the con- 
served quantities. By ergodicity the field config- 
urations sampled at large time intervals will be 
distributed according to the microcanonical en- 
semble. For large systems this is equivalent to 
the canonical ensemble: time average rs canonical 
average. Such averaging implies thermal equilib- 
rium but we can still study small departures from 
equilibrium by looking at linear response func- 
tions (O(i)O'(O)). Once the properties of these 
quantities are well understood we can turn to 
larger departures from equilibrium. 

In the canonical description in the temporal 
gauge A = 0, Gauss' law 5S/6A&, = = 
has the status of a constraint. The hamiltonian 
takes the form H = K + W, with now 
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and the canonical partition function 

Z= J DEDttDUD<p[Y[5(G*)} exp(-H/T). 
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By equipartition any initial configuration which is 
smooth on the lattice scale evolves into a 'rough' 
configuration, with temperature given by 
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where V is the volume and — 3 reflects the 
three Gauss contraints per lattice site x. The 
energy density diverges as the lattice spacing 
a — > at fixed temperature, which is the noto- 
rious Rayleigh-Einstein- Jeans divergence ||. 

3. Relation to dimensional reduction 

The REJ divergence is not the only one, which 
can be seen by studying time-independent cor- 
relation functions like (Q x (t)O y (t)) . It is then 
illuminating to integrate out the momenta by 
reintroducing Aq into the canonical partition 
function, using J exp^A^G^dA^ to represent 
<5(Gx)- This leads via a rescaling of A to the 
form of a three dimensional euclidean field the- 
ory, 
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Hence, for static quantities the classical theory 
is a dimensional reduction approximation to the 
quantum theory. 

Dimensional reduction is an accurate approxi- 
mation to a weakly coupled quantum theory in 
equillibrium at high temperature ||. The re- 
duced 3D theory is superrenormalizable which 
means that only mass counterterms are needed to 
obtain a finite theory. However, the above Sdr 
lacks a mass term for the 'adjoint Higgs field' Aq, 
since this is prohibited by the locality and gauge 
invariance of the classical action. In contrast, the 
usual derivation of the DR theory leads to addi- 
tional terms of the form 

Va TtA o and A^Tr^. 
The coefficient A^ turns out to be negligible, but 
the fact that [ia = in the above Sdr means that 



the there is no possibility for an A^ mass coun- 
terterm. It follows that the Debije screening mass 
mo associated with Aq is divergent in perturba- 
tion theory, mn oc 1/^/a. This need not be a 
disaster as a — > 0, because Aq will simply decou- 
ple as it gets heavy. Indeed, Ao is often integrated 
out explicitly as an additional approximation to 
dimensional reduction, because the renormalized 
Debije mass is large in perturbation theory. 

The classical action evidently has to be inter- 
preted as an effective action. Its parameters can 
be found by comparison with the dimensionally 
reduced quantum action with Aq integrated out, 
which is known analytically through perturbative 
calculations. This means that g 2 , A and /it 2 are 
explicitly known in terms of the corresponding 
parameters of the quantum theory, and a and T 
(|. In fact, g 2 w djjgi^)- This comparison also 
shows that the ratio ze/z^ » 1 to a very good 
approximation. 

We conclude, that for static observables the 
classical theory can approximate the quantum 
theory well. There is one undetermined parame- 
ter z = ze = z^. This parameter sets the time 
scale relative to the lattice spacing and it there- 
fore plays a role in dynamical quantities, to which 
we turn next. 

4. Dynamics 

Time dependent quantities like (O(t)O(0)) can 
be computed in the microcanonical ensemble by 
solving the equations of motion on a computer. 
We can 'help ergodicity' by using the canonical 
ensemble for generating many initial configura- 
tions and averaging over these. Luckily, there are 
now two good solutions to the algorithmic prob- 
lem of the implementation of the Gauss constraint 

To see how well the classical theory may fare 
for dynamical quantities, one may turn to pertu- 
bation theory. The problem of solving perturba- 
tively the equations of motion and averaging over 
initial conditions with the canonical ensemble has 
been studied in [QJ|] for 4 scalar field theory. 
It was concluded in Q that the mass countert- 
erm needed to make static correlation functions 
finite is also sufficient for obtaining finite time- 



3 




(a) 



(b) 



Figure 1. Gauge boson self energy diagrams. 



dependent correlation functions. The classical </> 4 
theory is renormalizable in this sense. Further- 
more, after matching parameters of the classi- 
cal theory to the dimensionally reduced quantum 
theory, the classical plasmon damping rate turned 
out to be identical to the quantum rate to lead- 
ing order in the temperature and coupling Q (see 
also J§]). Further analysis led to the conclusion 
that the classical A</> 4 theory can approximate the 
quantum theory for momenta and frequencies up 
to O(VAT), with corrections O(VA) @. 

The procedure of solving equations of motion 
and averaging over initial conditions is awkward 
analytically. A more convenient formalism can 
be given which is a classical analogue of thermal 
field theory [0. Alternatively one may use the 
imaginary time formalism of finite temperature 
quantum field theory and replace (after analytic 
continuation to real time) the Bose distribution 
by its high temperature approximation: 

1 1 T 

2 + e"/T-l w' ( J 

We assume this latter procedure for a brief dis- 
cussion of the situation in gauge theories. 
Consider the gauge boson selfenergy 

n M „(p,p°) = n$(p,P°) + n$(p,p°), as given 
by the diagrams in fig. [l], using the abelian Higgs 
model as a simple example. Both contributions 
(a) and (b) are linearly divergent (the quadratic 
divergence of the quantum theory is reduced to a 
linear divergence in the classical theory by (||)). 
Because of gauge invariance we may expect can- 
celations. Indeed, for external frequency p° = 
the fields are static and therefore the spatial com- 
ponents Ilmn(P) 0) and Ilmn(p, 0) are related by a 
Ward identity of the dimensionally reduced the- 
ory. As a consequence, their sum II m „(p, 0) is 



finite. However, for p° ^ 0, the cancellation is in- 
complete and H mn (p,p°) turns out to be linearly 
divergent, in both real and imaginary parts. This 
may seem unfamiliar because in the quantum the- 
ory at zero temperature the Ward indentity re- 
duces a quadratic divergence into a logarithmic 
divergence for all momenta. The reason such re- 
duction does not follow here is that at finite tem- 
perature U mn (p,p°) is not analytic at p = for 
p° 7^ 0. This is related to the physical process 
of Landau damping and a detailed physical pic- 
ture has been developed by Arnold which is 
also valid in the nonabelian case (on the lattice, 
measure-effects end up only into the momentum 
independent ITS). 

The momentum dependence of the divergent 
terms is complicated, typically of the hard ther- 
mal loop form ln[(p° + |p|)/(p°— |p|)]. This means 
that the divergencies in classical gauge theory are 
nonlocal in spacetime (in contrast to scalar field 
theory), which suggests nonlocal counterterms. 
This does not look attractive, although it may be 
possible to introduce new degrees of freedom to 
re-express such counterterms in local form. The 
situation is complicated, however, by the fact that 
on the lattice the divergencies lack rotational co- 
variance |^Jil|]. 

A different point of view prevails, in which the 
lattice spacing is supposed to stay finite, of order 
of the inverse temperature (the linear divergence 
oc g 2 T/a is cured to oc g 2 T 2 in the quantum the- 
ory), and where the hard thermal loop physics 
is to be added 'by hand'. There are difficulties 
with double counting, but practical proposals (l2| 
are being pursued. Yet another approach consists 
of deriving an effective theory in which high fre- 
quency modes are subdominant, such that regu- 
lator effects vanish as the lattice spacing goes to 
zero 1 13 . 

In conclusion, for time-dependent quantities 
the classical model suffers from hard thermal loop 
effects which have lattice artefacts and which are 
divergent as the lattice spacing goes to zero. This 
is also relevant for Lyapunov exponents (cf. the 
brief remarks in |I| ) . Awaiting a solution of these 
problems, we still may expect to obtain useful in- 
formation with the classical model at finite lattice 
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spacing, provided the results are interpreted with 
care 

5. Sphaleron rate 

The sphaleron rate, plays an important role in 
theories of baryogenesis It is the real time 
analogue of the topological susceptibility in imag- 
inary time, 
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in continuum notation. The usual expectations 
have been 

T/V 



aw = —, 
4ir 



oc exp[-£ s (T)/T], T<T C , 
w const, T > T c , 
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where E S (T) is the temperature dependent 
sphaleron energy (of order 10 TeV for T = 0, van- 
ishing near T c ) and const = 0(1). The prefactors 
of the exponential sphaleron suppression (||) have 
been calculated analytically and the main prob- 
lem for the numerical simulations is the compu- 
tation of k above T c . Recent work in pure SU(2) 
gauge theory Q and in SU(2)-Higgs theory j| 
using the same numerical implementation found 
k ~ 1.1 in the high temperature phase. Further- 
more, the results did not appear to depend on 
the lattice spacing. This is quite remarkable: a 
reduction of a by a factor 0.6 caused T/V to fall 
by factor (0.6) 4 « 0.13, in both phases J§. How- 
ever, the expected sphaleron suppression was not 
observed in the low temperature phase, and in- 
deed the results are considered to be wrong. The 
reason is that the imperfections of the 'naive' lat- 
tice implementation used for Tr FF imply not 
only the need for a multiplicative renormaliza- 
tion, but also a subtractive renormalization for 
the rate; these were not taken into account. Such 
a subtraction has been controversial in euclidean 
lattice QCD, but in practise it appears to work 
well || . 

Recently, there have been two new computa- 
tions using improved methods for obtaining T 
|l6|-|ll| , but before presenting these it is useful to 



mention new theoretical analysis questioning the 
lorical aw independence of 'const' in (Q). Arnold, 
Son and Yaffe (ASY) proposed the following pic- 
ture g§: 

the nonperturbative processes important for T 
occur on the typical momentum scale g 2 T; how- 
ever, due to Landau damping, the corresponding 
frequencies are not of the order g 2 T, but g 4 T. 

Hence, the prediction is T/V ss const' g 2 (g 2 T) 4 
or k, oc g 2 as g 2 — ► 0, in contrast to the lore that 
k in (Q) is practically independent of g 2 . 

Now g 2 T is the natural classical (dimensional 
reduction) scale for static processes (H/T = 
J d 3 xTr F mn F mn /&g 2 T H ), but the extra fac- 
tor of g 2 can only be supplied in a given regular- 
ization. On the lattice this involves the combina- 
tion (3 = 4/g 2 aT (recall (|l|)). So the prediction is 
that k oc i.e. oc a! Physically, the damping 
in the classical theory diverges as a — > 0, as dis- 
cussed in the previous section, and then nothing 
moves anymore. 

A detailed study of lattice effects in the clas- 
sical theory by Arnold jll| suggests that not all 
is lost: we should readjust the time scale in the 
classical theory proportional to the ratio of damp- 
ing rates of the quantum and classical theory. He 
obtained the following matching relation: 



k w 3.2 aw P K c i a 
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with an estimated error of about 30 % due to the 
cubic anisotropy of the lattice. 

The new methods for the computation of T use 
cooling (Ambj0rn and Krasnitz (AK) JI^pl ) or 
computing winding numbers of gauge transforma- 
tions (Moore and Turok (MT) Q). With either 
method the rate in the low temperature phase was 
found too small to be measureable, which is what 
one expects because of the sphaleron suppression 
@. Fig. H summarizes the results for K c i ass in the 
high temperature phase, in the form of /3re c i ass as 
a function of . The lower MT data are pertur- 
batively corrected for some discretization errors. 
The AK data should be compared with the upper 
MT data, and they can be seen to be compatible, 
but the trend as a function of (3 is different. If 
K c i ass indeed vanishes proportional to a we expect 
/S^ciass to approach a constant as (3~ l — > 0. This 
is clearly suggested by the MT data which can be 



25 



20 



15 



10 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 



Figure 2. Data for f3K c ] ass versus (3 1 . The 
straight lines are fits to the MT data fll6[ : 
Kciass = 25.1 /3 _1 - 56.7 /3 -2 (upper) and K c i ass = 
23.6 /T 1 - 78.1/3" 2 (lower). The AK data are 
connected by lines to guide the eye; lower: ref. 
0, upper: ref. @. 

fitted very well by a straight line (the slope repre- 
sents 0(a 2 ) corrections), but not by the AK data, 
for which K c iass itself and not /3K c iass is approxi- 
mately constant. Given the much larger leverage 
in the MT data one may tentatively conclude that 
Kciass ~ Of course, smaller /3 _1 are desirable 
to confirm this behavior, or not. However, small 
volume AK results also show K c i ass oc p7| . 
Using (g) and the corrected MT data the tenta- 
tive conclusion is then k « 76 %. 

6. Autocorrelation functions 

Apart from the correlation (|]) used for the 
computation of the sphaleron rate there is 
very little nonperturbative information on time- 
dependent correlation functions. Correlators 
Co = (O(t)O(O)} of gauge invariant fields 



O 



ip^ip = if, 



iTr^(d k -iA k )4>r a = W^, 



(<f> is the matrix version of the doublet ip) have 
been computed at zero momentum |2C[| . Their in- 
terpretation may be guided by perturbation the- 
ory. In the Higgs phase H and Wj? are equivalent 
to ip and Afr, to lowest order, but in the plasma 
(high temperature) phase the interpretation is 
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Figure 3. The autocorrelator Cij(t) in the Higgs 
phase for T/T c = 0.85, aT = 0.96 (small ampli- 
tude) and T/T c = 0.78, aT = 0.41 (large ampli- 
tude). The time is in lattice units. 



not so straightforward. In particular the anal- 
ogy with the confinement phase of QCD, which 
holds for time-independent correlation functions 
and which suggests pole dominance by confined 
states, appears to be misleading in this time- 
dependent case. Furthermore, there is to my 
knowledge no analogue of the zero temperature 
'theorem of the arbitrariness of the interpolating 
field'. 

Fig. U shows an example of Ch m the Higgs 
phase. Similar results were obtained for Cw(t)- 
From the damped oscillations 'plasmon' frequen- 
cies and damping rates lj and 7 were extracted 
by fitting to a pole dominance form C(t) = 
i?exp(— -ft) cos(ujt + a) + background. In the 
plasma phase the correlators are noisier; pole 
dominance is questionable and the extraction of 
7 is difficult especially for H. The resulting to 
and 7 appear to be approximately independent 
of the lattice spacing. Closer examination in the 
Higgs phase reveals, however, that the data for 
lu\v are compatible with a weak lattice spacing 
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ing wall between the Higgs and plasma phase, 
i] = pressure/ wall velocity, from the fluctuations 
of the wall position via a fluctation-dissipation ar- 
gument. Another computation was the change in 
J FF under influence of a chemical potential re- 
lated fermion production, in moving walls. There 
were several other interesting results and tech- 
niques, but there is no more space here for an 
illustration. 

I thank Gert Aarts, Peter Arnold and Alex 
Krasnitz for useful discussions. 



Figure 4. Zero momentum gauge field correlation 
function in Coulomb gauge, Ca(^)/Ca(0) versus 
t/0 2 , for /3 = 10, 12, 14; ref. @. 

dependence expected from the perturbative 1/y/a 
divergence. On the other hand, the damping rates 
are not divergent to leading order and both ju 
and "fw turn out to have magnitudes similar to 
the analytic results in the quantum theory 

Ambj0rn and Krasnitz | |l7j studied correla- 
tors of the gauge field itself in Coulomb gauge, 
O — > AjJ, for various momenta p, in pure SU(2) 
gauge theory. These are better suited for com- 
parison with perturbation theory and currently a 
hot issue is if the ASY analysis Jll],[l!| applies. 
This suggests overdamped correlation functions 
(which do not oscillate), while the classical diver- 
gence will show up in time scales increasing oc f3, 
i.e. in lattice units oc /3 2 . The analysis applies to a 
regime Pq < p 2 < T 2 , which is hard to satisfy in 
current simulations. Surprisingly the characteris- 
tics of overdamping and lattice time scale oc /3 2 
shows up even for p = 0, cf. Fig. |J. This remains 
to be explained. 

We have not shown here the richer structures 
found in the correlation functions in the plasma 
phase [ po|fl7j| , which are currently being analyzed 

nil- 

7. Dynamics of the phase transition 

In a very stimulating paper [p3| , Moore and 
Turok studied the real time properties of the elec- 
troweak phase transition by numerical simulation. 
They obtained the drag coefficient for the mov- 
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